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1 Introduction 



The universal properties of multiparton matrix elements in the infrared (soft and collinear) 
limit play a relevant role in our capability to make reliable QCD predictions for hard- 
scattering processes |], 0. 

At leading order in the QCD coupling as, these properties are embodied in process- 
independent factorization formulae of tree-level 0-|5[ and one-loop |6|-[8], [|] amplitudes that 
are at the basis of (at least) three important tools in perturbative QCD. The leading- 
logarithmic (LL) parton showers, which are implemented in Monte Carlo event generators 
IjH to describe the exclusive structure of hadronic final states, are based on these fac- 
torization formulae supplemented with 'jet calculus' techniques || and colour- coherence 
properties |]T0 , Analytical techniques to perform all-order resummation of logarithmi- 



cally enhanced contributions at next-to-leading logarithmic (NLL) accuracy |T^j rely on 
the factorization properties of soft and collinear emission. More recently, the leading-order 
factorization formulae have been fully exploited to set up completely general algorithms [7|, 



3~3| p~5[] to handle and cancel infrared singularities when combining tree-level and one-loop 
contributions in the evaluation of jet cross sections at the next-to-leading order (NLO) in 
perturbation theory. 

Any further (although challenging) theoretical improvement of these tools requires the 
understanding of infrared-factorization properties at the next order in as- 

This topic has recently received considerable attention. The infrared singular behaviour 
of two-loop QCD amplitudes has been discussed in Ref. JT^]. The collinear limit of one-loop 



amplitudes is known |T7[. In the case of tree- level amplitudes, the factorization structure 



in the double-soft [18| and triple-collinear [19| limits has been studied. 



In this letter, we reconsider the triple-collinear limit examined by Campbell and Glover 
19fl . In particular, we extend their results by fully taking into account azimuthal (spin) 



correlations. Besides improving the general understanding of the collinear behaviour of 
tree amplitudes, this extension is essential to apply some general methods to perform 
exact fixed-order calculations at the next-to-next-to-leading order (NNLO) in perturbation 
theory. For instance, the subtraction method [14, 15[] works by regularizing the infrared 



singularities of the tree-level matrix element by identifying and subtracting a proper local 
counterterm. Thus, the knowledge of the azimuthally averaged collinear limit studied in 



Ref. [|19j is not sufficient for this purpose. 



The outline of the paper is as follows. In Sect. |2], we define our notation and the collinear 
kinematics and, after reviewing the known collinear-factorization formulae at (D(as), we 
present their generalization at the next perturbative order. Our explicit results for the 
collinear-splitting functions at 0{a\) are given in Sect. |3|, where the comparison with 
the azimuthally-averaged case considered in Ref. [IS] is also discussed. Details on our 
calculation of the splitting functions are not discussed here and will be presented elsewhere 
f20|| . Some final remarks are left to Sect. |j. 
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2 Notation and kinematics 



We consider a generic scattering process involving final-state^ QCD partons (massless 
quarks and gluons) with momenta pi,p 2 ,... Non-QCD partons (7*, Z°, W , . . .), carry- 
ing a total momentum Q, are always understood. The corresponding tree-level matrix 
element is denoted by 

M^-f^^p,,...) , (1) 

where {ci,c 2 ,...}, {si,s 2 , •••} and {ai,a 2 , ■ ■ ■} are respectively colour, spin and flavour 
indices. The matrix element squared summed over final-state colours and spins will be 
denoted by \J^A. ai ,a 2 ,...ij ) iiV2, ■ ■ -)\ 2 - To the purpose of the present paper, it is also useful to 
consider the sum over the spins of all the final-state partons but one. For instance, if the 
sum over the spin polarizations of the parton a\ is not carried out, we define the following 
'spin-polarization tensor' 



spins ^s\,s\ colours 



(2) 



In the evaluation of the matrix element, we use dimensional regularization in d = 4 — 2e 
space-time dimensions and consider two helicity states for massless quarks and d — 2 helicity 
states for gluons. This defines the usual dimensional-regularization scheme. Thus, the 
fermion spin indices are s — 1,2 while to label the gluon spin it is convenient to use the 
corresponding Lorentz index [i = 1, . . . , d. The <i-dimensional average of the matrix element 
over the polarizations of a parton a is obtained by means of the factors 

T^ss' (3) 
for a fermion, and (the gauge terms are proportional either to p^ or to p u ) 

-JZ^ d ^P) = 2(1 _ e ) (~ 9 » v + gauge terms ) ( 4 ) 

with 

- Q^d^p) = d-2, p» d^ip) = d^(p) p u = 0, (5) 
for a gluon with on-shell momentum p. 

The relevant collinear limit at 0(as) is approached when the momenta of two partons, 
say pi and p 2 , become parallel. Usually, this limit is precisely defined as follows: 

p^zjf + k^- 1 ^^— , p£ = (1 _ - jfe£ _ -^1 — 



z 2p ■ n 1 — z 2p ■ n 

k 2 

S12 = 2pi ■ p 2 = 7- ± . k± — ► . (6) 

z(l - z) 

In Eq. (Q) the light-like (p 2 = 0) vector p^ denotes the collinear direction, while is 
an auxiliary light-like vector, which is necessary to specify the transverse component k± 

§The case of incoming partons can be recovered by simply crossing the parton flavours and momenta. 
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(k 2 ± < 0) (k ± ■ p = k A 



n 



0) or, equivalently, how the collinear direction is approached. 



In the small-fc^ limit (i.e. neglecting terms that are less singular than 1/Ar[), the square of 
the matrix element in Eq. ([!]) fulfils the following factorization formula [|l|] 



\M ai ,a 2 ,...(pi,P2,- ■ ■) 



Sl2 



4n^asT a s Up : 



(7) 



where /z is the dimensional- regularization scale. The spin-polarization tensor T^ 3 ' (p, . . .) is 
obtained by replacing the partons a\ and a 2 on the right-hand side of Eq. (|2]) with a single 
parton denoted by a. This parton carries the quantum numbers of the pair a>i + a 2 in the 
collinear limit. In other words, its momentum is p M and its other quantum numbers (flavour, 
colour) are obtained according to the following rule: anything + gluon gives anything and 
quark + antiquark gives gluon. 



The kernel P ai a 2 m Eq. ([!]) is the <i-dimensional Altarelli-Parisi splitting function [2I|. It 
depends not only on the momentum fraction z involved in the collinear splitting a — > ai+a 2 , 
but also on the transverse momentum k± and on the helicity of the parton a in the matrix 
element \'' s '"'(p, ■ ■ ■)■ More precisely, P aia2 is i n general a matrix acting on the spin 
indices s, s' of the parton a in the spin-polarization tensor T* s {jp, . . .). Because of these 



spin correlations, the spin-average square of the matrix element .Mq"' ;s ' 
simply factorized on the right-hand side of Eq. ([?]). 

The explicit expressions of P ai a 2 i f° r the splitting processes 

a{p) -> ai {zp + k ± + 0{k\)) + a 2 ((l - z)p -k ± + 0{k 2 ± )) , 

depend on the flavour of the partons ai, a 2 and are given by 

'1 + z 2 



P™(z,k ± ;e) = P™(z } k ± ;e)=5 ss ,C F 
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1 - Z 



e(l-z) 



cannot be 



(8) 



(9) 



P™'(z,k r ,e)=6 ss ,C F 



i + a-z) 



ez 



(10) 



P%(z,kr,e)=P£(z,kr,e) = T J 



P£(z,kr,e) = 2C A 



-^r +4z(i-^) fc±fc± 



A i 



-g 



Z 1 — Z 

+ 



,1 — Z z 
where the SU(N C ) QCD colour factors are 

TV 2 - 1 

C F ~- 



-2(l-e)z(l- z) 



hj I A; I 



C A = N C , T r = - 

2N r 2 



(11) 
(12) 

(13) 



and the spin indices of the parent parton a have been denoted by s, s' if a is a fermion and 
/i, v if a is a gluon. 

Note that when the parent parton is a fermion (cf. Eqs. flj^) and (|T0"|)) the splitting 
function is proportional to the unity matrix in the spin indices. Thus, in the factorization 
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formula flj]), spin correlations are effective only in the case of the collinear splitting of a 
gluon. Owing to the fc^-dependence of the gluon splitting functions in Eqs. QTTJ) and fll^D , 
these spin correlations produce a non-trivial azimuthal dependence with respect to the 
directions of the other momenta in the factorized matrix element. 

Equations (p|)-([T^) lead to the more familiar form of the <i-dimensional splitting func- 
tions only after average over the polarizations of the parton a. The rf-dimensional average 
is obtained by means of the factors in Eqs. @ and @). Denoting by (P ai a 2 ) the average of 
Pa x a 2 over the polarizations of the parent parton a, we have: 



P qg (z;e)) = (P qg (z;e)) =C 



1 + z 2 



e(l-z) 



(14) 



(P gq (z;e)) = (P gq (z;e)) =C h 



1 + 



€Z 



(15) 



(P qq -(z;e)) = {P- qq {z-e)) =T R 
(P gg (z;e)) = 2C A 



2z(l-z) 
1 - e 



z 1 - z . 
+ + z(l-z) 



1 - z 



(16) 
(17) 



In the rest of the paper we are interested in the collinear limit at 0(a\). In this case 
three parton momenta can simultaneously become parallel. Denoting these momenta by 
P11P2 and ps, their most general parametrization is 



Pi 



XiP^ + k" ±i - 



k 2 



Xi 2p ■ n 



i = 1,2,3, 



18) 



where, as in Eq. (g), the light-like vector p M denotes the collinear direction and the auxiliary 
light-like vector n M specifies how the collinear direction is approached (k±i-p = k±i -n — 0). 
Note that no other constraint (e.g. = 1 or k±i = 0) is imposed on the longitudinal 

and transverse variables X{ and k±i. Thus, we can easily consider any (asymmetric) collinear 
limit at once. 

In the triple-collinear limit, the matrix element squared \M. ai ,a2,a 3 ,...(pi,P2,P3, ■ ■ -)\ 2 has 
the singular behaviour |-A/f ai jCl2 a3j 

...(Pi>P2,P3, • • Ol 2 ~ V( ss ')> where s and s' can be either 
two-particle (sy = (pi + Pj) 2 ) or three-particle (si 23 = (pi + p 2 + P3) 2 ) sub-energies. More 
precisely, it can be shown [|l!|, ^(J that the matrix element squared still fulfils a factorization 
formula analogous to Eq. (|7|), namely 



\M ai ,a2,a 3 ,...(Pl,P2,P3, 



6 123 



(An^asY TZ(P, ■ ■ •) PZ 



0.20.3 



(19) 



Likewise in Eq. (|7|), the spin-polarization tensor T^ s ' (p, . . .) is obtained by replacing the 
partons ax, a 2 and 03 with a single parent parton, whose flavour a is determined (see Sect. |3|) 
by flavour conservation in the splitting process a — > a\ + a 2 + CI3. 
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The three-parton splitting functions P ai a 2 a A generalize the Altarelli-Parisi splitting func- 
tions in Eq. (|7|). The spin correlations produced by the collinear splitting are taken into 
account by the splitting functions in a universal way, i.e. independently of the specific 
matrix element on the right-hand side of Eq. (|T9|). Besides depending on the spin of the 
parent parton, the functions P ai a 2 a 3 depend on the momenta Pi,P2,P3- However, due to 
their invariance under longitudinal boosts along the collinear direction, the splitting func- 
tions can depend in a non-trivial way only on the sub-energy ratios Sij/s^ and on the 
following longitudinal and transverse variables: 



(20) 



K = ¥Li-^— J2 k 1j , (2i) 

l^k=l x k j = i 

which automatically satisfy the constraints Z)f =1 Zi — 1 and J2i=i K = 0. To simplify the 
explicit expressions of the splitting functions, we find it convenient to introduce also the 
variables 

t .. k = 2 Z ^ k ~ Z ^ k + Z -i^± Si . . (22) 

Zi ~\~ Zj Z'l ~\~ Zj 

The results of our calculation are presented in the next section. 



3 Collinear splitting functions at 0{al > 

To evaluate the three-parton splitting functions, we use power-counting arguments and the 
universal factorization properties of collinear singularities. The method |22] consists in 
directly computing process-independent Feynman subgraphs in a physical gauge. Details 
on the method and on our calculation are given in Ref. ||20|| . In the following we present 



the complete results for the spin-dependent splitting functions. 

The list of (non-vanishing) splitting processes that we have to consider is as follows: 

Q -" Qi + Q2 + Is , (q -* tfi + <?2 + ?s) , (23) 
q -> qi + q2 + «?3 , (q -> q\ + q-i + qz) , (24) 
q -»• 91 + 92 + q3 , (q -> 9\ + 92 + 93) , (25) 

9 -> 9i + 92 + q3 , (26) 
9 -> 91 + 92 + 93 ■ (27) 

The superscripts in q', q' denote fermions with different flavour with respect to q, q. The 
splitting functions for the processes in parenthesis in Eqs. ( p3|) and ( 24[) can be simply 



obtained by charge-conjugation invariance, i.e. Pq' iq ' 2 q 3 = Pq[q' 2 q 3 an d Pq iq2 q 3 — -Pgi? 2 ?3- ^ n 
summary, we have to compute five independent splitting functions. 

In the case of the splitting processes that involve a fermion as parent parton (see 



Eqs. (|23f)-(P5D), we find that spin correlations are absent. We can thus write 

3&* = 6SS ' ( P ^s) , (2? 
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and likewise for Pqf q2q3 and Pgf 92q3 - This feature is completely analogous to the 0(as) case 
and follows from helicity conservation in the the quark-gluon vector coupling. 

The spin-averaged splitting function for non-identical fermions in the final state is 



S123 

S12 



-12,3 , 4^3 + [Zl - Z 2 ) 



S12S123 



Z\ + z 2 



1 - 2e) [zi + z 2 



Sl2 

S123 



(29) 



The analogous splitting function in the case of final-state fermions with identical flavour 
can be written in terms of that in Eq. (p9"|), as follows 



where 



4^) = +(2^3) 



/p_(id) v 



(30) 



/p_(id) v 
X- 1 919293/ 



+ 



CV ( C*f — 



5123 I" 1 + z l __ 2z 2 
Sl2 [ 1 - Z 2 1- Z 3 

l + zl 



2s 



23 



g 123 Z l 
S12S13 2 



:i-^)(i-^s) 



S12 

1-30 

1-22 



+i+^i 



'1-Z 3 



2z 2 

1-^8. 



+ (2^3) . (31) 



The splitting function of the remaining quark-decay subprocess can be decomposed 
according to the different colour coefficients: 



(4,293) =C|(4 a 9 b 2 ) 93) + C F C A (P^l) 



919293' 

and the abelian and non-abelian contributions are 



(32) 



/ p(ab) \ 
\ 919293' 



+ 



S 2 

6 123 



2S13S23 



2.3 



l+i| 
ZlZ 2 



z\ + z\ 



Z\Z 2 



e(l + e) 



S123 



z 3 (l - Zi) + (1 - z 2 ) 



S13 

+ (1-0 



^2 

v l-e) — 

«13 



+ e 2 (l + z 3 ) - e(zl + z lZ2 + z\ 
(1^2) , 



22 



^2 



(33) 



/ p(nab) \ 
\ 919293/ 



+ 
+ 
+ 



^12,3 
M2 



1 e \ 

4 ~ 2 J + 2s 12 s 13 



? 23 r(l- Z3 )2(i_ e ) + 2z 3 



22 



z 2 2 (l- e ) + 2(l-z 2 ) 



1-^3 



'123 



4S1 3 S 2 3 



23 



:i-^3) 2 a 



2% 



2^2 



+ 6(l-e) 



g!23 
2S12 
■S123 

2s 13 



^ _ ^ zi{2 - 2zi + zj) -z 2 (Q- Qz 2 + zl) 2e z z{zi - 2z 2 ) - z 2 



'l-e) 



^2(1 - z 3 ) 



■\- Z2 f + zl-z 2 f 2(1 - z 2 )(z 2 - z 3 ) 



22(1 - z 3 ) 



— e 



z 2 (l ~ z 3 ) 



Ml - z 3 j 
- Zi + z 2 



Z\Z 2 



ZlZ 2 



+ (1^2) 



(34) 
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In the case of collinear decays of a gluon (see Eqs. (26, 27)), spin correlations are highly 
non-trivial. 



The colour-factor decomposition of the splitting function for the decay into a qq pair 
plus a gluon is 

A/ii/ _ /~i rp pM^(ab) n rp p /ii/(nab) /nr\ 

9i<?2?3 ~ L "F ± R r giQ2q3 ^ ^ R ^ gi q 2 q 3 , \<JO) 

where the abelian and non-abelian terms are given by 



P 



fiu (ah) 



919293 



+ 



-g 

43123 
S12S13 



_2 _|_ 2S123S23 + (1 — e)(si23 — S23)' 
S12S13 

^/l^/l2 H - k^k^ (1 c^jk^k 



V 1 

1 y ' 



(36) 



Pj/ii/ (nab) 
919293 



+ 



+ 



1 J s 123 

4 



6 23 



5 



23,1 



16- 



Z 2 Z 2 



^123 
S12S13 
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S123 Zi(l-Zi) 
2Sl23<r - 4(^^3 + W 
Sl23 1—^3 



Sl23 



- 2e) + 2 
-2si 23 S 



z 2 z 3 J \z 2 23 
- (1 " 

Si23 1 - 2Ti + 2z 2 



S12 2;i(l - Zl) 
., ^2(l-2^i) 



S23 2l(l - Zi) 



S12S23 

+ 4(^2^3 + ^3 ^2/ 



Zl(l - 2i) 
2^2(^3 - gi) 
Zl(l - ^1) 



+ (1 



2l(l - Z\, 

• + (2<->3) . 



8(1 - e)A#Jfe 



(37) 



In the case of gluon decay into three collinear gluons we find 
'fl-e 



919293 



r 2 



4s 2 



12 



-^2,3 + 16S123 



z\z\ 



1 - e)<T + 



S123 



4 



-5 



23 



2(1 



W - 2:3 
23 





-I)" 






\Z 2 






2 J . 



4zf 



S12S13 
S123 



:9 



i- uj 



S12 

4,2 2 2 3 + 2^l(l - Zi 



f ^3^3 



1-2^(1-^3) 
^l(l - Zi) 

-2z 2 \ 



1 



2:2(1 - z 2 ) 
1 - 2^(1 



2(1 



(l- Z2 )(l- Z3 ) 

'222(1-22) \ 



2223 / 

+ (5 permutations) . 



(3c' 



The splitting functions in Eqs. (|3 
the parent gluon according to Eq. 



2 3 (1 - 23) 

Wj can be averaged over the spin polarizations of 



(Pa 1 a 2 a 3 ) — 



2(1-6) 



0.20,3 



(39) 
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Performing the average we obtain 

,p(ab)_ v 
\ "'019293/ 



-2 - (1 - e)s 23 

$123 



Sl2 







+ 2 S?23 f 


\Sl2 


Sl3/ 





1 + 2z x + e - 2 



21 + z 2 



S123 



S13 



-, , 2 Z l + 2^2^3\ 

1+2 ' — r^J 

Z\ + 23 

1 + 2z 1 + e - 2-= 



(40) 



/ p(nab) v 



+ 



+ 



b 23,l 



+ 



'123 



^123 

2s 13 

2^ 

1 e 
4 + 2 



22 



2Si 3 S 2 3 
1 + 



23 



2z 3 (1 - 23 - 2ziz 2 ) 



2l(l 



Zl) (1-6)^(1-^0 
222(1-22) 



1 + zf 
zi(l - z l) 



+ 



Zl (l- Zl ) (1 - £ )z 1 (l - Zl ) 

21(23 - 2 2 ) 2 - 2z 2 z 3 (l + Zi) 



(1 - e) Zl (l - z x ) 



'123 



2Si 2 Si 3 



( 1 + , ; _£l±i^)} + (2 „3) 



(41) 



\P9i92g3) 



+ 



+ 




4<? 2 

^ 6 12 



"12,3 



+ t(1-6) 



(l-^3(l-2 3 )) 2 
Z 3 Zi(l - 2l) 

1 + 2^(1 + ^) 
2(l-z 2 )(l-z 3 ) 



+ 



'123 



_ Sl23 ^ z l z 2 ~ 1 _|_ 2lZ2 " 
Sl2 [ 1 — 2 3 Z 3 

ziz 2 (l - 2; 2 )(1 - 2z 3 ) 



2 3 5 
- + 2 + 2 Z3 



S12S13 
1 - 2zi(l 



23(1 - 23) 



+ z 2 z 3 - 2 + 



2z 2 z 3 



+ (5 permutations) 



zi(l + 2zi) 



(42) 



The (9(ag)-collinear behaviour of tree-level QCD matrix elements has been indepen- 



dently examined by Campbell and Glover ||19|| . Their study differs in many respects from 



our analysis. Taken for granted the universal factorization formula (|T9|), they compute the 
three-parton splitting functions by directly performing the collinear limit of the explicit ex- 
pressions of the 7* — ► four- and five-parton squared matrix elements. Moreover, they treat 
the colour structure in a different way and consider the collinear limit of the colour-ordered 
sub-amplitudes 0. Finally, they neglect spin correlations and present only the explicit 
expressions of the polarization-averaged splitting functions. 

By properly taking into account the differences in the colour treatment, we have com- 
pared our results with those of Ref. and found complete agreement^ for the spin- 
averaged splitting functions. To be precise, the colour-connected splitting functions P ai a 2 a 3 ^c 
of Ref. are related to our spin-averaged splitting functions as follows: 



( P <?i 9^,93 ) 
/p_(id) v 

/ p(ab) \ 
\ 919293' 



' 123 ^j^f^Jj-, pnon— ident. 



939^2 ->« 



c F c 



'123 

4 

s 2 

'123 p 



•> ident. 



— C a I P 1 ' - 

2 / 939192-^9 ' 



^Note that the published version o f Ref. fl9|| contai ns two misprints (in Eqs. (5.8) and (5.19)) that have 
been corrected in the archive version hep-ph/9710255 v3. 
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s 2 1 

/ p(nab) \ fl23 ^ / p , p - P - - ) (A$) 

\ 913293/ 4 2 ^ 939192->g ^ r q:i929i^q r q:i9i92^q) > l^ d J 

„2 

/p(ab) v ^123 

\ r giq2qsl ^ r q2~giq^9 > 

2 1 

; p(nab) v _ 5 123 _ ( p _ _|_P_ _P_ "| 

\- r 9i9293/ ^ 2 v 919392— >9 ~T ^939291^9 ^929193^9/ ' 

(^919293) = [^919293-9 + ( 5 permutations)] . 

Owing to the completely different methods used by the two groups, this agreement can be 
regarded as an important cross-check of the calculations. 



4 Summary 

We have considered the three-parton collinear limit of tree-level QCD amplitudes. In 
this limit the singular behaviour of the matrix element squared is given by the universal 
factorization formula fll9D and is controlled by process-independent splitting functions, 
which are analogous to the Altarelli-Parisi splitting functions. In Sect. |3| we have presented 
the explicit expressions of the splitting functions at 0(a|), taking fully into account spin 
correlations. 

These splitting functions are one of the necessary ingredients needed to extend QCD 
predictions at higher perturbative orders. In particular, they are relevant to perform an- 
alytic resummed calculations beyond NLL accuracy and to set up general methods to 
compute jet cross sections at NNLO. The knowledge of the collinear splitting functions, 
when combined with a consistent analysis of soft-gluon coherence, can also give prospects 
of improving the logarithmic accuracy of parton showers available at present for Monte 
Carlo event generators. 
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